eo CHINESE JOURNAL ree 
hag OF ENGINEERING MATHEMATICS asad 


Article ID:1005-3085(2010)01-0161-07 


On the Solution of a Sylvester Equation 
AV + BW = EVF Appearing in Singular 
Control Theory* 


GONG Wen-zhen 
(Department of Mathematics and Computer Science, Yulin Normal University, Yulin 537000) 


Abstract: This paper deals with the solution of the matrix equation AV + BW = EV F associated to 
a linear singular systems and subject to some rank and regional pole-placement constraints. 
Under the [-strongly controllability of the singular system, a sequence of coordinate trans- 
formations is proposed such that the considered problem can be solved through a Sylvester 
equation associated to a controllable reduced-order normal system. The solution of the 
Sylvester equation can also be found by the regional pole-placement in LMI regions of the 
complex plane. The results provide great convenience to the computation and analysis of 
the solutions to this class of equations. 
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1 = Introduction 


Consider the singular system (FẸ, A, B} described by the following equation 


Et = Az + Bu, 
(1) 
y = Cr. 


Here z € R”, u € R”, E, AER", BE R”*?, C e R™™*", rankE =q, rankB = p < q and 
rankC = m. As in singular systems (1), many problems, such as observer design!™?], and the 
problem of eigenstructure assignmentl3-®}, are closely related with the matrix equation of the 
following form 

AV + BW = EVF. (2) 


Here A, B, E are given matrices, while F € R(I-P)*(9-P); V e C™*(-P) and W e CP*(9-p) 
are to be determined; F is in the Jordan form with arbitrary given eigenvalues. In [7], Eugenio 
and Vilemar have obtained an explicit solutions for (2). Duan!*~4] gave a complete analytical 
restriction free parametric solution of (2), and then presented a new solution of (2). These 
existing solutions are applicable in some problems like eigenstructure assignment. However, 
many strict constraint conditions are imposed and large computational load are needed to 
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obtain these solutions of (2) in the above literature. In this note, a new and simple approach 
to solve matrix equation (2) is developed. As a by-product, we show that this matrix equation 
(2) corresponds to the solution of some known Lyapunov like equations that allows to obtain 
regional pole-placement in LMI-type regions of the complex plane. The structure of this paper 
is as follows. The next section describes some background material and specific mathematical 
tools needed for the algorithm development. Section 3 presents the proposed solution of matrix 
equation (2). A necessary and sufficient condition is first derived for the existence of the solution 
for (2). Then an algorithm is presented to compute the solution of (2) by the eigenstructure 
assignment. In section 4, a numerical example is presented. In section 5, a conclusion ends the 


paper. 


2 Preliminaries 


Assume that matrix E € R”*” is singular and rankE = q < n, then there exist two 
nonsingular matrices P, Q with appropriate dimensions such that 
I, 0 Ai. A B 
Eo =PEQ=| * des PAG = 7P |, Bo=PB=| ' 
0 0 Ag, Ao Bo 
Let I’ be a region in the open left half complex plane, T C C~, symmetric with respect to 
the real axis. Let L € R"*("-® is any full column rank matrix satisfying EL = 0. Let 


L 
oiee aaa 
L2 
; ee Lı 
according to this partition, EL = 0 becomes Ec = 0, which leads to Lı = 0 and 


2 
L € R©-9)x@-4) is an arbitrary nonsingular matrix. 
Definition 1 System (1) is said to be T-strongly controllable if and only if the following 


conditions are satisfied 
rank(AE - A B)=n, VAEC, AgT, (3) 


rank(E AL B)=n. (4) 


3 Main results 


3.1 The solution of matrix equation 
Let 


B, = (A12 Bı), B= (Az Bə), 

ie B R V; 

R= E € Re* atp), V = PEV = ; , T=(0 Taa); 
Bo 0 
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we obtain the following Theorem. 

Theorem 1 Assume that the matrix B is of full column rank, that is, rankB = n — q+D, 
the necessary and sufficient condition for the existence of the solution to the matrix equation 
(2) is that system (1) is I-strongly controllable. 

The proof of theorem1 needs the following Lemma 1. 

Lemma 1 The system (1) is I’-strongly controllable if and only if 


AI-—Ac B 
rank =2n—q, VAEC, AGT and (Ag. Bo) =n-—™m. 
T 0 


Proof To establish the sufficient and necessary conditions, we consider the following equa- 
tions 
AIT—Ai —Ai2 Bı 
rank(AEZ — A B)=rank(\Ec — Ac Bo) = rank (5) 
-A2 Áz Bə 
and 
ei Ig ÁrLl2 Bı 
rank(E AL B)= rank(PEQ PAQQ™'L PB) = rank . (6) 
0 AzL2 Bz 


Definition 1, (5) and (6) imply Lemma 1. 

Proof of Theorem 1: The proof is constructive. 

Sufficiency If system (1) is T-strongly controllable. Under the assumption that the matrix 
B is of full column rank, we can choose a matrix R € R”*(9-P) such that rank(R B) =n. Let 


z M: 
(R B= 
Nı 
From Lemma 1, we have that 
rank(Ag2 B2) = rank, = rankTĒ = rank(TR TB) = rankT =n — q. (7) 
Let 
X = M,AcR-—M,AcB(TB)*TR, (8) 
Z = Mı = AcB(I - (TB)*(TB)), (9) 
= Mı 0 
R BO I 0 0 
U, = , U=] a a z ~ |o Us=|N, OG}, (10) 
0 0 I (TB)+TR I-(TB)+(TB) (TB)t ea 


where At denotes any generalized inverse of matrix, satisfying AA+A = A. Consider the 


following equations 
| M-Ac B | M—MAcR M,AcB 
rank 


TR TB 


n — q + rank 
T 0 


ll 


2n — 2q + p + rank(AI —- X Z). 
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From Lemma 1 and (8), we get 
rank(I-—-X Z)=q-p, VAEC, AGT. (11) 


From (9), for any given F € R(&@-?)*(9-P), o(F) CT, we can choose the appropriate matrix Y, 
such that 


F=X+ZY. (12) 
Taking 
K = (TB)'TR - (I ~ (TB)+(TB))Y, (13) 
then we have TBK = TR, or 
T(R- BK) =0. (14) 
(13) is written as 
F = M, AcR -— M,AcBK. (15) 
Let 
Kı = MAcR-— N, AcBK, (16) 


By combining (13), (16) and (17), we get 


RF + BK, + AcBK = AcR, (17) 
or 
Ac(R — BK) -(R- BK)F = BJ, (18) 
where 
Ji zako = 
J=- =K,+KF, J e RO-9*@-?), Jy e RP*OP), (19) 
2 
= Ry _(n— 
(R—-BK)= | | R; € RTO, (20) 
0 


Substituting (17) and (18) into (16) yields 
Ay Ry Aj2J1 + By Jo Ry 
+ = F. (21) 
Ag Ry Az22J1 + B2J2 0 


V=Q >, WH=d. 


1 


Let 


From (20), we get that 
PAQQ71V + PBW = PEQQ™'VF, (22) 


which prove the sufficiency of Theorem 1. 
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Necessity From the proof of sufficiency, we get that F = X + ZY CT if and only if (9) 
holds, and the fact (9) holds if and only if 


M-Ac B 
rank =2n-q, VAEC, AGT and (Ag. Be) =n-gQ, 
T 0 


and combining Lemma 1, we easily obtain the proof of necessity. 

3.2 Algorithm 

Let the conditions of Theorem 1 be satisfied, the algorithm for the solution of the matrix 
equation (2) can be described as follows. 


I 0 
Step 1 Compute matrices P, Q such that PEQ = | , then compute 
0 0 


A A B 
Ac = PAQ = u 42 |p =PB= All OVE 
Aoi A22 Bo V2 


Step 2 Choose a matrix R such that rank(R B) =n, compute (R B)-! and let 


(R By = | 
Nı 
Compute X and Z from (8) and (9), respectively. For given F C I, using the pole placement 
method in the context state space systems to compute the gain matrix Y, such that F = X+ZY, 
then compute K, and K from (14) and (17). 
Step 3 Compute 


J z 
J=- | | =K,+KF, J € RCX), Jye RPXG-?), Ry = (Iq 0)(R- BK). 


H 
V=Q », Wshh. 
Ji 

3.3 Regional pole placement in LMI regions 

In this section, we will discusse the solution of matrix equation AV + BW = EVF, such 
that o(F) CO = {A € C | fe = A+ AL + AL7}, where A = AT e RO&?)*(G-P) and 
L= LT e R-P)» (a-p), 

Theorem 2 Assume that system (1) is T-strongly controllable, then there exist a solution 
to AV + BW = EVF with o(F) € © if and only if there exist a symmetric positive matrix 
M e RG&-»)x(@-P) and a matrix N € R&~?)*(@-?) such that 


A®M+L@(MXT)+L@(NZ7)+L7 8(MXT) +L" (NZT) <0, (23) 
where 


X = M,AcçR— M,AcB(TB)*TR, Z= Mı = M,AcB(I - (TB)*(TB)), 


166 CHINESE JOURNAL OF ENGINEERING MATHEMATICS VOL. 27 


and & is the Kronecker’s product. In this case F = X + ZY with Y = N7M~-T. 
Proof It holds that o(F) C © if and only if o(F7) C ©. Since o(F7) C O, there exists a 
symmetric positive definite matrix M € R(1-P)x(4-P) such that (see [8]) 


A®M+L@(MF")+L7(MF") <0, (24) 
From the above result, we have, and by substituting this value into (25), and by putting 
N = MYT , we obtain the desired result. 
4 An illustrative example 


Consider the matrices Æ, A, B and C defined below, corresponding to a I’-strongly-controllable, 
whose finite poles are given by o(E, A) = {j, j}. 


100 0 0 01 0 
0100 0 -1 0 0 

E= , A= , B= , B= 
0010 -1 1 0 0 1 
0000 1 0 0 1 -1 


Let P = Q = I, then Ec = E, Ac = A, T = (0 0 01). From (8) and (9), we obtain 


[to] [2] 


For any given F CT, for example, let 


Compute 
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Let 

0 -i1 

R 

V=Q| “|= , W=h=(-1 1). 

Jy 0 1 

-1 2 
ne f -1 0 ; . 
Now it is easy to illustrate that, for F = , V, W is the solution to the matrix 
0 - 


equation AV + BW = EVF. 


5 Conclusion 


In this paper, the solution of the matrix equation AV + BW = EVF associated to singular 
systems has been settled by a new method. A numerical example has shown that the new 
method can provide an interesting framework to implicitly and optimized solutions for the 
considered constrained Sylvester equation. This solution of this problem can also be obtained 
by regional pole placement in LMI regions. 
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KFT MABE HA Sylevster H42 AV + BW = EV F WARR 
48M He 
(EAT ZARA SABLE AR, “PK 537000) 

HB: ELZETTE MARA Sylevster JFE AV + BW = EVF RRM, AEEA SHAK 
FE REAR RRR AC AS. ZED N RAE TSR REE BIN PAU ER, BEAT 
38 52 AY Sylevster Fy #2 AS AE AY fal RUE Ly ES) R BE HE AEN Sylevster F7 FEAR. Sb, PT 
Æ FEY Sylevster FERA th Bt ARAI LMI RRR ECR RAH. BCE RMT A 
PIXE A FEY RS RAB A 7 E. 

Ki: MAB: Sylevster FH: KARA 


